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Precise Distances for Main-Belt Asteroids in Only Two Nights 

Aren N. Heinze^’^ and Stanimir Metchev^’^’^ 

ABSTRACT 

We present a method for calculating precise distances to asteroids using only two 
nights of data from a single location — far too little for an orbit — by exploiting the 
angular reflex motion of the asteroids due to Earth’s axial rotation. We refer to this 
as the rotational reflex velocity method. While the concept is simple and well-known, 
it has not been previously exploited for surveys of main-belt asteroids. We offer a 
mathematical development, estimates of the errors of the approximation, and a demon¬ 
stration using a sample of 197 asteroids observed for two nights with a small, 0.9-meter 
telescope. This demonstration used digital tracking to enhance detection sensitivity 
for faint asteroids, but our distance determination works with any detection method. 
Forty-eight asteroids in our sample had known orbits prior to our observations, and for 
these we demonstrate a mean fractional error of only 1.6% between the distances we 
calculate and those given in ephemerides from the Minor Planet Center. In contrast 
to our two-night results, distance determination by fitting approximate orbits requires 
observations spanning 7-10 nights. Once an asteroid’s distance is known, its absolute 
magnitude and size (given a statistically-estimated albedo) may immediately be calcu¬ 
lated. Our method will therefore greatly enhance the efficiency with which 4-meter and 
larger telescopes can probe the size distribution of small (e.g. 100 meter) main belt 
asteroids. This distribution remains poorly known, yet encodes information about the 
collisional evolution of the asteroid belt — and hence the history of the Solar System. 

Subject headings: astrometry, celestial mechanics, ephemerides, minor planets, aster¬ 
oids: general 
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1. Introduction 


The main asteroid belt is a relic from the formation of the Solar System. Although much of its 
mass has been lost, it retains a great deal of information about Solar System history and presents 
us with a laboratory in which we can study collisional processes that once operated throughout the 
circumsolar disk in which Earth and the other planets were formed. One of the most straightforward 
observables constraining such processes is the asteroid belt’s size-frequency distribution (SFD; 
Bottke et al. 2005a). The current mai n belt’s SFD can b e successfully modeled as the result of 
4.5 billion years of collisional evolution (jBottke et al.ll2005al : Ide Elia fc Brunini 1120071 ). While such 
models fit the ‘collisional wave’ set up by 100 km asteroids able to survive unshattered through the 
age of the Solar System, they cannot be observationally tested in the 100 meter size range. 


Objects in this size range are very interesting, because they supply most near-Earth asteroids 
and meteorites by shattering one another and/ or migrating inward via Yarkovsky and resonance ef¬ 
fects ( Farinella et al .11 1994 : Bottke et al.ll2005b l. Modern 8-10 meter telescopes can detect them, but 
monitoring them over many nights to determine an orbit requires a prohibitively large time invest¬ 
ment for such powerful telescopes (e.g., 7-10 nights; Gladman et al. 2009). Thus their distances 


and sizes remain unknown, and detailed an alyses are confined to larger objects ( Gladm an et al 


2009) or use only rough statistical distances ( Yoshida et al. 2003 : Yoshida Sz Nakamura 2007l i. 


We present a method to obtain precise distances to main belt asteroids (MBAs) using only 
two nights of observations. Distances translate directly into absolute magnitudes and hence to 
sizes given a reasonable assumption for the albedo distribution. This method, which we refer to 
as rotational reflex velocity (RRV), will greatly increase the efficiency of surveys aimed at probing 
collisional evolution in the Solar System by measuring the SFDs for extremely small MBAs. 



^The WIYN Observatory is a joint facility of the University of Wisconsin-Madison, Indiana University, the National 
Optical Astronomy Observatory and the University of Missouri. 
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2. Distances from Rotational Reflex Velocity 

Suppose that at a given instant, an asteroid located a distance d from an Earth-based observer 
is moving with velocity Va, while the observer is moving with velocity Vq (e.g., the orbital velocity 
of the Earth). The angular velocity at which the observer sees the asteroid move relative to distant 
stars is given by: 


Va± - Vo± 

where the T subscript indicates the vector component perpendicular to the line of sight, so that 
Va± — Vo± is the projection of the asteroid’s relative velocity onto the plane of the sky. 

Although Vo can be accurately calculated for any Earth-based observation, the velocity Va of a 
newly discovered asteroid is always unknown initially, and hence the distance cannot be calculated 
by simply plugging the measured value of uj into Equation [TJ Given appropriate measurements, 
however, we can isolate the component of uj that reflects the observer’s motion around the geocenter 
due to Earth’s rotation, and from this calculate the distance. This is the essence of the RRV method 
for distance determination. 

The velocity Vq of an observer on the surface of the Earth can be expressed as the sum of 
Earth’s orbital velocity Vorb and the velocity Vrot with which the Earth’s axial rotation carries the 
observer around the geocenter. Neglecting the slight asphericity of the Earth, Vrot = Veq cos 9, 
where 6 is the observer’s terrestrial latitude and Vgq is the Earth’s equatorial rotation velocity of 
1674.4 km/hr. For convenience, we define Vrei as the asteroid’s velocity relative to the geocenter: 
Vrei = Va — Vorb- The angular velocity ujg that would be measured by an ideal observer located 
at the geocenter then depends only on Vrei and the distance, but the angular velocity Uo that is 
measured by a real observer based on the Earth’s surface depends also on Vrot- The two angular 
velocities are given by: 


( 1 ) 


UJ 


9 ~ 


VrelA- 

d 


( 2 ) 


Vrell. Vrotl. 

UJo = - 1 - 

d 

If we could measure ujg, we could therefore calculate the distance: 


( 3 ) 


, Vrot 

d = - 

UJg - UJo 

where we have dropped the T subscript, because it will henceforward apply to all physical velocities 
in our calculations. 


( 4 ) 
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Now suppose that the asteroid is observed near midnight on two different nights, that the two 
observations are separated by exactly one sidereal day, and that the position and angular velocity 
Wo are recorded for each observation. The angular distance the asteroid moved between the two 
observations will thus be accurately known; call this (j). Because exactly one full rotation of the 
Earth elapsed between the two observations, the observer’s position relative to the geocenter is the 
same for both of them. Thus, the average geocentric angular velocity of the asteroid in between 
the two measurements is uog = 4>/At, where At is the elapsed time between the observations: one 
sidereal day. 

Let the measured values of Wo on the first and second nights be Woi and LO 02 , and similarly 
let the perpendicular rotational velocities (which are obtained by calculation, not measurement) 
be Vroti and Vrot 2 - We can then evaluate the difference between geocentric and observer angular 
velocities twice: the average of Ug — Woi and Ug — uJo 2 will be a factor of \/2 more precise than a 
single measurement if the uncertainty on ujg measurement is much smaller than on Woi and Wo 2 - 
This is likely to be the case, since ujg is based on a longer temporal baseline. The distance is then 
given by: 


, Vrotl T Vrot2 

“ = 7 -^-7-7 (a) 

(Wg - tool) + (Wg - UJ02) 

So far we have assumed that Vrei and d show no appreciable change over the 24 hour period 
of the measurements, so that Wg — Woi and Wg — Wo 2 are effectively two measurements of the same 
quantity. We will now determine the errors that result when this assumption is violated. We will 
first consider changes in Vrei and d that are linear in time: that is, when the first time derivatives 
Vrei and d are significant but the second derivatives Vrei and d are not. 

We will parameterize the change in Vrei by Au = {vrei 2 — Vreii)l‘^, and the change in d by 
e = {d 2 — di)/{2d). For convenience, we will also use <5 = {vrot 2 — Vroti)/{‘^Vrot) since the observations 
do not have to be taken exactly one sidereal day apart, and hence the Earth’s projected rotational 
velocity may vary. For well-optimized observations <5 should always be close to zero. Note that 
while e and 6 are unit-less, fractional changes, Av still has units of linear velocity. We then obtain: 


Wg — UJol 


Vrei _ At?—(1 —(5)tJrot 

d (l—e)d 


A'g UJo2 


Vrei _ Vrel+^V-{'^+S)Vrot 

d (l-|-e)d 


( 6 ) 


We substitute these into Equation [SJ Retaining terms only to second order in the variational 
quantities Au, e, and 5, after a good deal of algebra we find that the fractional error on the 
calculated distance is: 


dcalc dtrue 2 ^rel Av o 

- - -= e-e- Veo-e 

^true '^rot '^rot 


( 7 ) 
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For main belt asteroids e is always very small (i.e. < 10“^). This renders 6 innocuous even 
though it can reach 0.3 for poorly-planned observations (Section 13.21) . The velocity change Av 
approaches zero for asteroids at opposition, where the accuracy of Equation [5] can be better than 
10~®, but even months from opposition the error remains only of order 10“^. Equation [5] is therefore 
very accurate under all conditions where the required input data could reasonably be obtained, and 
the uncertainties even of highly accurate measurements (Section 13.3j) will generally be larger than 
the equation’s intrinsic error. 

We caution the reader against a much less accurate alternative to Equation [5] that may appear 
more attractive because it is obtained simply by applying Equation [J] individually to each night 
and then averaging the resulting two distance measurements: 

d = ( 8 ) 

— ^ol — Wo2 

Despite its intuitive simplicity, this equation should never be used. Its fractional error expansion 
contains the term which has very large values for objects more than a few days from 

opposition, as illustrated by EigurelH By contrast, for Equation [5] the precursors of the Av^ 
term and several other error terms cancel in the denominator, leaving only the very much smaller 
error terms given by Equation [71 Equation [5] should therefore be avoided, and Equation |S] (or its 
generalized form. Equation [T8|) should always be used for RRV distances. 


3. Asteroid Angular Velocities from Positions Measured at Discrete Times 
3.1. Basic Equations and Methodology 

Eor simplicity, the discussion in Section [2] assumed the observer could make instantaneous 
angular velocity measurements. In practice, of course, angular velocity measurements must be 
made over a period of time. This has very little effect on our mathematical derivation. The only 
change is that the projected velocities of the observer relative to the geocenter, Vroti and Vrot 2 , 
should not be considered as instantaneous but rather as averages over the same period of time in 
which the angular velocities lOoi and u}o 2 were measured. 

Positions and angular velocities of asteroids are measured in terms of right ascension (RA) and 
declination (DEC) on the celestial sphere, and thus are two dimensional. Eor example, an observer 
could measure an asteroid as having an angular velocity of -37.36 arcsec/hr in RA (the negative sign 
indicating westward rather than eastward motion) and 10.32 arcsec/hr in DEC, with the positive 
sign indicating the object is moving toward celestial north. In principle, independent versions of 
Equation [5] could be constructed for each dimension, with angular velocities in RA corresponding 
exactly to projected east-west motion of the observer and DEC velocities corresponding to projected 
north-south motion. However, the true rotational velocity of the observer is always strictly east- 
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Days from opposition 



Fig. 1.— Errors of RRV distances calculated from JPL Horizons ephemerides for representative 
known MBAs. Black points are for Equation 1 181 a version of Equation [5] that is generalized for non- 
optimal timing of observations (Section 18.2|) . while gray points are for a similarly generalized form 
of Equation [8l Left: Near optimally timed observations with x (Equation [T7]) equal to -0.0007, 
corresponding to t 2 m — iim equal to one minute less than a sidereal day. Right: Non-optimally 
timed observations with y = 0.052, corresponding to t 2 m — tim equal to 1.3 hours more than a 
sidereal day. Equation [18] yields accurate distances even with poorly timed observations, while 
Equation [8] should never be used. ‘Noise’ in the plots is due to roundoff error in the ephemerides 
rather than intrinsic error in the equations. It is more severe in the right-hand case because of the 
shorter temporal baseline on each night (1.7 hr vs. 2.9 hr). 
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west, and although its projection Vrot can have a north-south component, this component is usually 
small — i.e., less than 10% of the east-west component for observations with mean hour angle in 
the range ±1.0 targeting asteroids within 20° of the celestial equator. As we will find below, for 
our test observations only the RA velocity components are large enough to supply useful distances. 
We expect this will generally be the case. 

In most asteroid surveys, many helds are observed, and each held is visited at least three times 
per night. Each visit yields a single celestial position of each detected asteroid. Two observations 
of each asteroid on each night will suffice for our purposes. Using the traditional symbols a and 5 
for RA and DEC, suppose that on night 1 we obtain observations aia,<5ia at time tia and aih,6ih 
at time tih, and similarly on night 2. Then we can also dehne mean positions and the mean time 
aim = (aia ± aib)l2, 5im = ((^la ± <5 i 6)/2, and Um = (tia ± tib)l2, with the analogous quantities 
being calculated for night 2. Note that the times here are measured on a continuous sequence and 
not reset when the date changes. Eor simplicity, at present we will suppose further that the time of 
each measurement on night 2 is exactly one sidereal day later than the time of the corresponding 
measurement on night 1: thus t 2 a—tia = t 2 b—tib = one sidereal day. The required angular velocities 
are then given by the following equations, where for simplicity we show only the RA component: 


^Ib ^la 
‘^ola — , , 

^Ib ^la 

(9) 

_ a2b — a2a 

^o2a. — , , 

^2b ^2a 

(10) 

a2m aim 

^2m ^Im 

(11) 


where the denominator of the last equation is of course also equal to one sidereal day. 

The only other quantities needed for the distance calculation are Vroti and Vrot 2 : the observer’s 
projected rotational velocities relative to the geocenter. As illustrated by EigureO the observer’s 
projected motion relative to the geocenter traces out a segment of an ellipse whose axis ratio is 
the sine of the asteroid’s declination. The projected distance of the observer from the geocenter is 
given by: 


x{t) = R© cos(6'zat) sin 

y{t) = -R© cos{eiat) sin((5) cos sm{9iat) 


( 12 ) 


where x and y are positive toward the east and north; Oiat is the observer’s latitude on the Earth; S 
is the declination of the asteroid; t is the time of the observation; is the time when the asteroid 
stands at zero hour angle for the observer; and ts is the length of a sidereal day. Note that the 
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Fig. 2.— Example illustration of an observer’s motion due to the Earth’s rotation, as viewed from 
an asteroid. Observation A and Observation B are both made from the same location on Earth’s 
surface (at 20° N, roughly the latitude of Mauna Kea), but they are separated by an elapsed 
time of 3.989 hr, corresponding to 60° of sidereal rotation. The rotational velocity of Earth’s 
surface at 20° N lat is 1573.4 km/hr. The effect on the asteroid’s measured positions is determined 
by the projected distance moved between the observations, which has eastward and northward 
components as labeled. The projected rotational velocity Vrot that should be used in the RRV 
distance calculation (Equation [5] or Equation I18p is therefore 5625.7 km/3.989 hr = 1410.3 km/hr 
eastward and 423.8 km/3.989 hr = 106.2 km/hr northward. 
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argument 27r(t — tz)/Ts is the hour angle at which the observer sees the asteroid: at t = tz the 
asteroid crosses the meridian and stands closest to the observer’s zenith and in the best sky position 
for observing. 

For observations at tia and tu, the projected linear velocity Vroti in the east-west (x) and 
north-south (y) directions is: 


_ - x{tia) 

'^rotlx — 

^Ib ^la 


(13) 


'^rotly 


yjhb) - yjtia) 

tlb tla 


(14) 


Figure [2] illustrates a specific example of two observations of an asteroid at -12° DEC, made from 
a single observing site at 20° N lat on the Earth’s surface. The first observation is made when the 
asteroid is rising, and the second about four hours later, shortly after it has crossed the meridian. 
The rotational velocity of Earth at 20° N lat is 1573.4 km/hr and is of course strictly eastward, but 
projecting the velocity into a plane perpendicular to the line of sight reduces the absolute magnitude 
and introduces a north-south component (except for asteroids exactly at 0° DEC). In this example 
the projected velocities east and north are 1410.3 km/hr and 106.2 km/hr. Since celestial RA and 
DEC correspond exactly to terrestrial longitude and latitude, the east-west physical velocity of 
Equation [13] can be combined with the RA components of the angular velocities from Equations 
MM and plugged directly into Equation [5] to yield the distance. 

We will now provide a simplified but basically realistic example. Suppose the MBA illustrated 
by Eigure[2]is being observed at opposition. It will be in its retrograde loop, and hence its motion 
in RA will be negative (westward). The angular velocity Wo measured by the observer will be faster 
to the west than the geocentric angular velocity Ug, because on the night side of the Earth the 
rotational velocity adds to the orbital velocity, making the asteroid fall back westward at a faster 
rate. Suppose Uga = —35 arcsec/hr and cooia = —36 arcsec/hr. The rotational reflex velocity, given 
by the difference ujga — ojoia, is then 1 arcsec/hr or 4.85 x 10“® rad/hr. Taking the physical and 
angular velocities to be exactly the same on the second night for simplicity. Equation [5] reduces to: 


d = 


'^rotx 


1410.3 km/hr 
4.85 X 10~® rad/hr 


= 2.91 X 10^ km = 1.94 AU 


(15) 


With Equation [15] we have calculat ed the distance to an asteroid based on the known prop¬ 
erties of Earth’s rotation. Interestingly, iBernstein fc Kushalani (]2000]) have developed a closely 
analogous method for determining the distance to a Kuiper Belt object using the known properties 
of Eart h’s orbit. Just as our d e rivati on does not (explicitly) include the effects of Earth’s orbit, 
that of [Bernstein &: Kushalani ([20001 1 does not include the effects of Earth’s rotation. Neglecting 
Earth’s orbit works in our case because, as we have demonstrated in Section [J] the dominant errors 
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due to orbital acceleration of both t he Earth and the asteroid can cel when the correct form of the 
distance equation is used. Similarly, iBernstein &: Kushalani ((20001) can safely neglect Earth’s rota¬ 
tion because the Kuiper Belt objects whose distances they aim to measure are sufficiently fa r away 
that their RRV signature is negligible. It follows that the method of lBernstein &: Kushalani (j2000|) 
should not be directly applied to measuring asteroid distances without some form of correction for 
the RRV signal, although (as they point out) it is useful for distinguishing true Kuiper Belt objects 
from asteroids that have similar angular velocities for a few days near their turnaround points. 


3.2. The Case of Observations that are Not Optimally Timed 

Up to now we have assumed for simplicity that the observations obtained on successive nights 
are exactly one sidereal day apart. While it is not difficult to plan and execute observations that 
satisfy this criterion to within a few minutes, clouds or instrument failures could intervene, or one 
could be processing an archival data set that was not taken with RRV distances in mind. In such 
cases the mean times tim and t 2 m aiay be separated by up to a couple of hours less or more than one 
sidereal day. A direct calculation of the geocentric angular velocity uig then becomes impossible. 
In its place, we must use instead the mean angular velocity across the two nights, defined by: 


^fia — 


(^2m 

^2m ilm 


(16) 


We will define three additional quantities. The time tg is the moment exactly one sidereal day 
before t 2 m (in the ideal case, tim would have been equal to tg)- The projected rotational velocity 
Vrot,s is the observer’s mean projected rotational velocity over the time interval between tim and 
tg. Finally, the fractional parameter y is given by: 


X = 


is ilm 
i2m ilm 


Note that y goes to zero in the case of optimally-timed measurements. 
The distance formula analogous to Equation [5] becomes: 


(17) 


, '^rotl '^rot2 ‘^X^rot,s 

d = 7-r-7-^ (18) 

{U)g — LOoi) -|- [LOg — U!o2 ) 

The error terms for Equation [18] are the same as those given in Equation [71 with the addition 
of an e^y term. For main belt asteroids, this will always be negligible. In closing our discussion 
of error in the case of non-optimal observation timing, we note that 6, the fractional change in 
projected rotational velocity from night 1 to night 2, is likely to be considerably larger in this 
case, up to 0.3. However, in the error terms for Equation [T8| it is always multiplied by e, which is 
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always of order 10 ^ or smaller for MBAs, and thus should not constitute a major source of error. 
Equation [18] is therefore the central formula for asteroid distance determination with RRV. 

In addition to our analytical error calculation (Equation [7|) , we have probed the errors of our 
approximations numerically using ephemerides of known objects from the JPL Horizons ephemeris 
generator. Note that no actual observations are involved in this test, only ephemeris positions 
and velocities for objects with well-known orbits. In Eigure [1] we plot the error of Equation [T8| 
as determined by this test. In contrast to our analytical error calculations, this test intrinsically 
includes the contributions of second and higher-order time derivatives of d and VreU since it uses 
JPL Horizons ephemerides for real objects. The plotted errors are affected by roundoff error in 
the ephemerides (accurate only to 0.15 arcsec in RA), which introduces pseudo-noise such that 
only the lower envelope of the plotted points can be meaningfully compared with our analytical 
error estimates. Nevertheless, the plots show that Equation [T8| produces accurate distances out to 
at least 60 days from opposition; that accurate distances can be obtained from observations that 
are not optimally timed; and that the second-order time derivatives Vrei and d do not introduce 
significant error. 


3.3. Measurement Precision Required for Good Distances 

Figured] although based on calculated ephemerides and not actual data, shows that relatively 
small position errors can create distance inaccuracies at the level of a few percent. We now consider 
the measurement precision that is necessary to yield distances to various levels of accuracy. 

In Equation 1151 the difference between the night-to-night mean angular velocity in RA (w^o) 
and the angular velocity in RA observed within a given night {oJoa) was one arcsec/hr, and this 
rotational reflex velocity corresponded to a distance of about 2 AU. This is typical for observations 
obtained near opposition at terrestrial latitudes of 20-30°. The fractional uncertainty of the distance 
is equal to that of — ^oa- This will generally be dominated by the uncertainty in Uoa, since 
it is calculated within a single night and is thus based on a smaller temporal baseline than 
Thus, if uioa is measured with an accuracy of 0.1 arcsec/hr, we should expect the distance to have 
an accuracy of only 10% at 2 AU. At 1 AU, the RRV signature is twice as large: uo^a — ^oa will 
therefore be about 2 arcsec/hr, and the calculation will be accurate to 5%. 

As we will now illustrate, one can usually obtain angular velocity measurements with better 
than 0.1 arcsec/hr uncertainties in practice. In our current paradigm, cjoa is based on two measure¬ 
ments of celestial positions separated by a time tf, — ta (i.e.. Equation or fTOl) . If the uncertainty 
on RA for these measurements is cTo, the uncertainty on ojoa will be: 

(Tq-V^ 

(Jtjj - 




(19) 
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Based on our experience, well-sampled (e.g. < 0.5 arcsec/pixel) images taken in ~ 1.5 arc- 
second seeing can yield aa = 0.03-0.05 arcsec for bright objects (e.g., those at least several times 
brighter than the lOcr detection threshold). For faint objects near the detection limit, uncertainties 
of cJq ~ 0.1 arcsec are more typicaS A temporal baseline of 5 hr is easy to obtain near opposition. 
Thus, we can expect = 0.008-0.014 arcsec/hr for bright objects and = 0.03 arcsec/hr for 
fainter objects. At 2 AU, distances will therefore be accurate to about 3% for faint objects and 
1-1.5% for bright objects, while at a distance of 1 AU, the accuracy should be about 1.5% even for 
faint objects. These accuracies are sufficient for good analyses of the size statistics of small MBAs, 
especially given that their albedos are unknown and must be tre ated as a statistical distribution 
spanning a factor of ~10 (jPravec et al.ll2012l : iMasiero et al.ll2013l ). Where the objects are faint or 


distant (e.g. 3-5 AU), or where extremely accurate distances are desired, it may be possible to 
obtain smaller values of iTq, by obtaining each RA measurement not from a single image but from 
the average over a set of images acquired close together in time. Extending the same principle 
further, we can use the technique of digital tracking to obtain accurate measurements for objects 
too faint even to be detected in individual images, as we describe below. 


4. Measurements Using Digital Tracking Data 
4.1. A Brief Introduction to Digital Tracking 


The observations we use to demonstrate the usefulness of rotational reflex distances for as¬ 


teroids were obtained using digital tracking, as described in our companion paper iHeinze et al 


(I 2 OI 5 I I. These data consist of 126 two-minute exposures of a single starfield that we acquired using 
the WIYN 0.9-meter telescope at Kitt Peak on the night of April 19, and 130 identically acquired 
images from the following night. 

Very briefly, a digital tracking search involves shifting and stacking such sets of images to reveal 
moving objects too faint to be detected on any individual frame. A separate trial stack is produced 
for each angular velocity vector in a hnely-sampled grid that spans the full range of possible sky 
velocit ies for the target population. For MBAs, a digital tracking data set can span only one 
night (jHeinze et al.ll2015l l. hence we analyze our April 19 and April 20 observations independently. 
Digital tracking has been used to great effect for Kuiper Belt objects (e.g. Allen et al. 2001, Fraser 
& Kavelaars 2009, and many others), but has not typically been used for faster-moving asteroids 
where a larger numbers of trial vectors must be probed. It is now computationally tractable even for 
asteroids, however, and produces a factor of ~ 10 increase in sensitivity over conventional methods. 
This sensitivity increase enables us to detect 215 asteroids within the 1-degree field of our test 
observations, despite using only a small, 0.9-meter telescope. We obtained precise angular velocity 


^Note that the distance uncertainty must be evaluated differently for digital tracking observations, where the 
angular velocity is directly measured (see Section 14.211 rather than being obtained from two positions. 
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measurements on both nights for 197 of these asteroids, including 48 previously known objects with 
accurate orbits. These last allow us to test the accuracy of RRV distance measurements using real 
data. 


4.2. Angular Velocity Measurement with Digital Tracking 

The geocentric angular velocity ojg of an asteroid detected on two subsequent nights is deter¬ 
mined in exactly the same way with digital tracking observations as with conventional data. The 
determination of the observed angular velocities oJq on each individual night is different, however 
— particularly because digital tracking enables the accurate measurement of asteroids that cannot 
even be detected on an individual image. 

Any asteroid detection in a digital tracking search occurs on a particular trial image stack, 
which corresponds to a particular angular velocity. Thus, for example, our automated digital 
tracking routine might search for objects with angular velocities between -50 and -20 arcsec/hr in 
RA and -10 and -|-20 arcsec/hr in DEC, using a grid spacing of 0.2 arcsec/hr, and might detect a 
particular asteroid as a bright point source on a trial stack whose shifts correspond to an angular 
velocity of -41.2 arcsec/hr in RA and 11.8 arcsec/hr in DEC. 

We improve on the relatively crude angular velocity measurement that is implicit in such 
a detection by probing a new, much more finely spaced grid of angular velocity vectors, using 
small postage-stamp images centered on the detected asteroid to make the search computationally 
tractable. For each trial stack in this finely sampled grid, we calculate the measured flux of the 
asteroid within a small aperture of radius roughly equal to the half-width of the point spread 
function (PSF) of our image. We then perform a 2-D quadratic fit to the measured flux as a 
function of the trial angular velocity. The angular velocity at which the quadratic ht reaches its 
peak value then constitutes an accurate measurement of Wq for that asteroid, and the uncertainty 
on ojo is derived from the uncertainty of the quadratic fit, which is naturally larger for faint objects 
that are more affected by sky background noise. 

All the inputs required to calculate the distance using Equation [18] are now available except 
Vrot- As in the case of discrete observations, Equation [12] indicates how Vrot should be calculated, 
but there is one subtle aspect. The appropriate value of Vrot is not a projected distance divided by 
an elapsed time as it was in the discrete case (Figure [2]), but neither is it the time average of the 
projected velocities x{t) and y{t). Instead, the value of Vrot that corresponds to the angular velocity 
measured by digital tracking is given by the slopes of the best linear fits to x and y as functions of 
time. Figure [3] illustrates the difference between linear fit velocities and average velocities. 

For our data, the linear fit velocities of the Kitt Peak observer were 1340.0 km/hr eastward 
and -4.8 km/hr northward on April 19; and 1337.3 km/hr eastward and -16.3 km/hr northward on 
April 20; by contrast the average velocities were 1286.6 km/hr eastward and -6.9 km/hr northward 
on April 19 and 1283.0 km/hr eastward and -14.6 km/hr northward on April 20. In this case, the 
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north-south velocities are negligible for practical purposes. 


5. Results from our Test Data 


5.1. Initial Distance Determinations 


Our observations on 2013 April 19 span a temporal range of 5.72 hr. Those obtained on April 
20 have the same range, and are centered 23.968 hr later in time. As this interval is only 2 minutes 
longer than a sidereal day, the timing of our observations is almost exactly optimal as defined in 
Section EU the parameter y from Equation [T7] is only 0.0014. We detected a total of 215 asteroids 
through digital tracking analysis of our observations in a field of view only slightly larger than 1 
square degree. We reported positions and brightnesses for all of these objects to the Minor Planet 
Center. Of these 215 MBAs, 202 were detected on both April 19 and April 20, and 197 of these two- 
night objects had sufficiently accurate measurements on both nights for good distance calculations. 
Of these, 48 were previously discovered objects with accurately known orbits, nine were previously 
discovered objects with poorly known orbits prior to our observations, and the remai ning 140 were 
new d iscoveries for which we received designations from the Minor Planet Center. 

(|2015l i provides further details. 


Heinze et al 


We have used Equation [18] to calculate distances for all 197 objects with sufficient measure¬ 
ments. The 48 objects in this sample with well-known orbits allow us to test the accuracy of 
rotational reflex distances using real data. The mean absolute error of our calculated distances for 
these objects is 0.042 AU, and the mean absolute fractional error is 2.2%. As expected, errors are 
smaller for more nearby objects. Distances with 2% accuracy are more than sufficient to analyze 
the size statistics of MBAs. Nevertheless, even better results are possible, as we will see below. 


5.2. Correcting for Bias from Track Curvature 

Close investigation reveals some evidence of systematic error in the distance calculations de¬ 
scribed in Section 15.11 The weighted average signed fractional error for the 48 known objects is 
—0.0179 ±0.0017: i.e., the calculated distances are 1.8% too small on average, and the offset is lOcr 
signihcant. 

It might be suspected that this bias is due to using the wrong value for Vrot — he., that we 
were wrong to use linear fit velocities rather than time-averaged velocities (Section 14.2p . This is 
easily demonstrated not to be the case, however. Substituting averaged values in place of linear 
fit values for Vrot changes the result too much and in the wrong direction: the mean systematic 
offset goes from -1.8% to -5.2%. The true cause of the bias is more subtle, and points to a further 
interesting measurement that can be extracted from our data. 
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Fig. 3.— The projected motion of a ground-based observer due to the Earth’s rotation, which affects 
the measured positions and angular velocities of detected asteroids. Left: Schematic diagram of 
Earth viewed from the perspective of asteroids in our April 19, 2013 field, which was near -12° 
declination. The heavy black curve is composed of 126 X’s marking the projected position (x, y 
from Equation I12p of the Kitt Peak (31.9583° N lat) observer at the moment each of our images 
was taken. Top Right: x{t) and y{t) for our April 19 images, showing the linear fit velocities 
that should be used for comparison to linear, constant-velocity digital tracking results. Bottom 
Right: The velocities x{t) and y{t) for our data, showing how the mean velocities are significantly 
different from the best-fit velocities shown in the top panel. 
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With digital tracking, our measurement of the single-night observed angular velocity oJq is 
made over a period of time equal to the temporal span of our observations on each night. The 
instantaneous angular velocity over this interval is continuously changing, with the rotational reflex 
velocity itself being the dominant source of change. For example, the rotation of the Earth causes 
the westward sky motions of the asteroids to be slightly faster near midnight than they are at 
the beginning and end of each night’s observations (see Figured]). As we stated in Section 
and illustrated by Figure E] we can account for the acceleration and curvature of an asteroid’s 
observed motion by using a value for Vrot that is obtained by a linear fit to the observer’s projected 
coordinates x{t) and y{t). 

This assumes, however, that the angular velocity measurement is effectively an unweighted fit 
to angular position as a function of time. Such is not necessarily the case. For example, the images 
taken near midnight, when the westward sky motion is at its fastest, may be more sensitive due 
to better seeing or lower atmospheric extinction. The digital tracking stack will then effectively 
weight them more highly and produce a mean angular velocity that is systematically too fast. This 
effect would produce a systematic underestimation of the asteroids’ geocentric distances: exactly 
what we observe for known asteroids in our data. 

The effective weighting of different images will not necessarily be the same for different as¬ 
teroids. An asteroid that passes too close to a bright star near midnight will only be measured 
near the beginning and end of the night, resulting in a measured angular velocity slower than the 
true average value. Asteroids can also exhibit significant brightness changes during the night due 
to their own rotation, which will change the effective weighting of different images in the final 
measurement of their angular velocities. 

If asteroids actually moved in straight lines at constant velocity, none of these weighting effects 
would introduce error into our angular velocity determinations. If we could predict the curvature 
and acceleration of each asteroid’s track due to the Earth’s rotation, and apply small shifts to each 
image in our digital tracking stacks in order to linearize the motion, we would therefore remove the 
bias in our angular velocity determinations. We now describe how to do this. 

The angular distances between the position an observer measures for an asteroid and the 
position it would have if measured from the geocenter are simply a reflection of Equation [T^ scaled 
by the distance: 


Aa{t) 


x{t) 

d 


A{(() = -VA 


( 20 ) 


Averaged over the span of a digital tracking integration, Aa{t) and A6{t) contribute a constant 
offset to the measured position of an asteroid, and their first time derivatives Aa{t) and A5{t) 
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produce a constant offset (i.e. the RRV signal itself) to its measured angular velocity. Neither of 
these constant offsets concerns us here: we wish to isolate and remove only the component that 
represents curvature and acceleration. Thus, we create new functions ^(t) and that behave 
exactly like Aa{t) and A5{t) except that over the span of our digital tracking integrations, their 
means and the means of their first time derivatives are zero: 


^{t) = Aa{t) - {Aa{t)){t - (t)) - {Aa{t)) 
C{t) = A6{t)-{ASmt-{t))-{A6{t)) 


( 21 ) 


where the notation () denotes the time-average of the enclosed quantity over the span of a digital 
tracking integration. The left panel of Figure S] presents Equation [21] in graphical form, plotting 
the curved track of an asteroid in ^ space over the course of our April 19 observations. The 
calculation of ^ and C of course requires knowledge of the asteroid’s distance d, but the 2% accurate 
values from Section 15.11 are more than sufficient for this purpose. 

Thus, we adjust the methodology described in Section 14.21 for measuring precise angular ve¬ 
locities from digital tracking data: we apply additional shifts equal to —^(t), —C(^) to each postage 
stamp image before the final stack. This straightens the asteroid’s track so that curvature can no 
longer influence the calculated angular velocity. The distance can be recalculated based on the 
resulting new, de-biased angular velocities. Note that because Equation [2T] subtracts average reflex 
velocities, we must now use averaged rather than linear-fit values for Vrot (Section 14.21 and Figure 
[3| illustrate how the velocities differ). 

Compared to the calculation in Section [5Tl our curvature correction reduces the mean absolute 
error of the distances for known objects from 0.042 AU to 0.032 AU, and the mean absolute 
fractional error from 2.2% to 1.6%. More significantly, the new calculation changes the weighted 
average signed fractional error from —0.0179 ± 0.0017 to —0.0028 ± 0.0017. The systematic offset 
is thus reduced by a factor of six, and the new value of -0.3% is statistically consistent with zero. 

The left panel of Figure [5] illustrates the precision of our measured distances for known objects. 
We emphasize that these distance measurements were entirely calculated from first principles and 
known quantities describing the Earth. We have used the known asteroids in our data to quantify 
our errors, but not to determine any of the parameters used in the distance calculation. 


5.3. Measurements of Track Curvature within a Single Night 

Beyond removing systematic errors in our measurements of asteroids’ mean angular velocities, 
we can in principle use Equation [2T] to measure the curvatures of asteroid tracks within just a single 
night’s data. This is an extremely challenging measurement because the curvatures are so small: 
the RMS curvature amplitude, defined as the root mean of C^(t) -|- C^(^)) is less than 0.1 arcsec 
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Fig. 4.— Left: The curvature terms ^(t),C(^) from Equation [2T] for an asteroid at a distance of 2 
AU in our April 19 data. The projected rotational velocity of the observer peaks at zero hour angle 
(near the center of the observing sequence); hence the asteroid’s retrograde (westward) angular 
velocity peaks at that time. Right: Measured rms curvature amplitudes (that is, the root mean of 
^^(t) + C^(t)) for asteroids brighter than R = 20.5 in our April 19 data, plotted against distance and 
compared to the predicted (not fit) relationship from Equations [20] and [2TJ Despite large scatter 
and a small systematic offset, the predicted 1/d dependence of curvature amplitude on distance is 
clearly seen, and the RMS error from this relationship is only 0.03 arcsec. In principle, curvature 
measurements with digital tracking data enable a measurement of an asteroid’s geocentric distance 
based on only one night’s data. 
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Fig. 5.— Left: The precision of our final RRV distance measurements for known objects. The line 
is the 1:1 correspondence and is not a ht: our distances are calculated from basic geometry and 
known quantities describing the Earth, not based on a fit to known objects in our data. As expected 
given the 1/d dependence of the RRV signal, our measurements are most precise for the nearest 
objects. Right: Histograms of the apparent magnitudes (and corresponding physical diameters) 
for all asteroids (solid line) and for previously known objects (dashed line) measured in our data. 
The current census of the main belt becomes substantially incomplete at a diameter of about 2 km. 
By contrast, we have detected dozens of new asteroids in the 200-500 meter size range. 
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for an asteroid at 2 AU. We attempt to measure it using our curvature-corrected postage-stamp 
analysis as described in Section [521 but instead of calculating ^{t) and ({t) based on a previously- 
calculated approximate distance, we probe a range of RMS curvature amplitudes for each asteroid. 
For each trial curvature amplitude we calculate ^{t) and Ci't), apply the appropriate shifts to the 
postage stamp images, and perform a quadratic fit to flux as a function of angular velocity as 
described in Section 1421 We compare the peak asteroid flux values produced by the quadratic fits 
for the different trial curvature amplitudes, and identify the curvature amplitude that yields the 
highest peak asteroid flux. This should be the curvature value that resulted in all of the individual 
asteroid images being most accurately registered. As the curvature amplitudes map uniquely to 
distances (Equations [20| and [2T]), our identification of the best-fit amplitude constitutes a distance 
measurement, in principle. However, the uncertainty of curvature amplitudes measured from the 
current data set is large enough that we do not choose to calculate individual curvature-based 
distances. 

Instead, we plot the best-fit curvature amplitudes as a function of the known RRV distances 
calculated in Section 15.21 A clear inverse relationship emerges, as predicted by Equation |20l The 
relationship gets tighter when we reject the minority of asteroids with out-of-range best-fit curvature 
values and all asteroids fainter than R = 20.5 mag, which cannot be measured as accurately. The 
48 objects meeting these criteria in our April 19 observations are plotted in the right-hand panel 
of Figure 01 which shows the RMS curvature amplitude as a function of geocentric distance in AU, 
and compares the data with the predicted curve from Equations [20] and EH Although there is a 
mild systematic offset, the slope of the best linear fit to curvature as a function of d~^ for these 48 
asteroids is 9cr significant, matches the predicted value to within 0.4a, and has an intercept only 
1.4(7 away from the predicted value of zero. The RMS error of the measured curvature amplitudes 
relative to their predicted values is only 0.03 arcsec. Thus we have effectively performed astrometry 
of moving, 20th magnitude objects with a precision of 30 milli-arcsec using only an 0.9m telescope. 

In principle, our success at measuring asteroid curvatures within a single night’s data means 
that digital tracking can be used to measure the geocentric distances to asteroids based on obser¬ 
vations from only one night. However, our current curvature measurements have uncertainties that 
are too large to yield useful distances for individual objects: we have been able to demonstrate the 
efficacy of our curvature-measuring methodology only by recourse to an ensemble of 48 asteroids 
with distances known by other means. With a larger telescope, more accurate curvature measure¬ 
ments of fainter asteroids would be possible, especially in good seeing. These could conceivably 
yield meaningful distance measurements for unknown main-belt asteroids over just a single night, 
albeit only for objects considerably brighter than the digital tracking detection limit. By contrast, 
the two-night RRV calculations of Sections 15.11 15.21 yield precise distances even for faint objects 
just above the detection threshold. There is therefore little reason to use the curvature method for 
asteroids measured on more than one night. 

However, where only one night’s data is available, distances based on digital-tracking curvature 
measurements could be very valuable. This is particularly true in at least two cases that may arise 
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in future surveys. The first is the case of an NEO survey where some objects — perhaps due to 
their fast motion — have inadvertently been measured on only one night. Being much closer to 
the Earth than MBAs, NEOs will show much larger and more easily measurable curvature, which 
is likely to yield useful distances. Such distances could enable the inclusion in a statistical analysis 
of objects that were otherwise unusable due to their one-night status. The second case is the one- 
night detection of a nearby NEO moving at nearly the same space velocity as Earth, such that 
its slow angular velocity mimics that of a much more distant object. This scenario is statistically 
rare but troubling, because such an object would likely be overlooked in an NEO survey and yet 
could be an incoming Earth-impactor. If the discovery survey used digital tracking, a curvature 
analysis would reveal the object’s actual, very small geocentric distance. We note in closing that 
curvature measurements can in principle be performed on asteroids detected with conventional 
methods rather than digital tracking, but they will normally be less accurate because of the smaller 
number of images available. 


5.4. Distances and Sizes of Detected Asteroids 


While our current data only suggest the possibility of curvature-based distance measurements 
from a single night, they allow precise two-night distance measurements from rotational reflex 
velocity (Section 15.21) . These precise distance measurements enable us to calculate the absolute 
magnitudes of each of our detected asteroids, and hence their physical diameters modulo the un¬ 
certainty in albedo. The histograms of these values are shown in the right-hand panel of Figure [5j 
At least half of the asteroids we have measured are smaller than 1 km. By contrast, only a small 
fraction of known MBAs are in this size range. The smallest objects we have detected are 300 
meters in diameter under the assumption of low, 5% albedo, but could be as small as 150 meters if 


their albedo is 25%. The statistics of MBAs in this size range are already known (jGladman et al 


20091 ) based on a sample several times larger than we present herein. Given this, the detailed com¬ 
pleteness analysis that would be required to convert our detections into a mea surement of the SEP 
is not worthwhile. The aims of the current work and our companion paper iHeinze et al.l (l2015l ) 
are instead the validation of the RRV technique and of digital tracking for asteroids, respectively. 
Similar observations using a larger telescope, however, would extend to much smaller objects and 
accurately measure the SFD of a previously unexplored size regime in the main belt. 


6. Conclusion 

We have described how the reflex angular velocity of asteroids due to Earth’s rotation can be 
used to determine the distances to main belt asteroids based on only two nights of observations. We 
refer to this as the rotational reflex velocity (RRV) method for measuring asteroid distances (Equa¬ 
tion [18]). The required approximations are accurate to about 10“^ (Equation [7|) , and measurement 
uncertainties are typically 1-3%. Such distances can be used to calculate precise size statistics of 
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small main belt asteroids using a much smaller investment of time on a large telescope than has 
previously been required. Accurate RRV distances may be calculated either from c onventional 


asteroid-search observat ions, from data analyzed by the technique of digital tracking (|Zhai et al 


2014 : iHeinze et al.ll2015l b or from any other specialized asteroid-observation technique (e.g. Milani 


et al. 1996 or Gural et al. 2005) that accurately measures the celestial coordinates of the asteroids. 


We have tested RRV distance determination with a data set acquired on April 19 and 20, 2013 
using the WIYN 0.9-meter telescope at Kitt Peak. Using measurements based on digital tracking, 
we have calculated distances to 197 asteroids in this data set. While the majority of these were new 
discoveries in our data, 48 of them are previously known objects with accurate orbits. These allowed 
us to test the accuracy of rotational reflex distances. A preliminary analysis yielded distances with 
a mean fractional error of only 2.2%. Without digital tracking, our detections would be confined 
to substantially brighter asteroids (or, the observations would require a larger telescope), but the 
accuracy of RRV distances for detected objects would likely be at least this good (see Section 1331) . 

We have identified a 1.8% systematic error in our preliminary RRV distances. We have linked 
this error to the curvature of the asteroids’ observed tracks that results from Earth’s rotation. 
This curvature, combined with the non-uniform sensitivity of images obtained under conditions of 
different seeing and airmass, causes a slight bias in the angular velocity measurements obtained from 
our digital tracking image-stacks. Finding our initial distance measurements more than sufficient 
to predict curvature amplitudes, we have calculated the form and amplitude of the curvature for 
each of our asteroids. We have applied a curvature correction to our image stacks to linearize the 
asteroid motions. Re-calculating the distances using angular velocities from curvature-corrected 
image stacks reduces the systematic bias from 1.8% to 0.3 ± 0.2%, and thus eliminates it as a 
significant effect. The corrected distances have a mean fractional error of only 1.6% for the 48 
known objects in our data. Note that such a curvature correction would not be needed for RRV 
distances based on asteroids detected on discrete images. The calculation of accurate RRV distances 
would therefore be simpler for conventional asteroid detections, although digital tracking can detect 
much fainter objects. 


Besides correcting our digital tracking angular velocity measurements for the curvature due to 
Earth’s rotation, we have attempted to measure this curvature by creating several different image 
stacks for each asteroid at a range of different curvature amplitudes. As the curvature amplitude 
is proportional to l/d, this constitutes a rotational reflex velocity measurement of an asteroid’s 
distance based on just a single night’s data. We find that we can indeed measure the curvature 
amplitudes for asteroids brighter than R = 20.5 mag with a precision of about 30 milli-arcsec. 
The measured curvature values for 48 asteroids brighter than this limit show the expected 1/d 
dependence at 9a significance, and the best-fit slope is within OAa of the predicted value. While 
distance measurements based on single-night curvature amplitudes are too noisy to be useful for 
asteroids in our current data set (and are not needed given our highly precise two-night values), 
single-night distances could be valuable in specific cases for future surveys. In particular, curvature 
measurements can identify nearby objects moving at slow angular velocities characteristic of a much 
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more distant population. 

Our precise distances from two-night RRV measurements allow us to calculate absolute mag¬ 
nitudes and hence approximate diameters for our newly discovered asteroids. Our faintest objects 
have Hpi ~ 21.5 mag and hence diameters of 130-300 meters depending on their albedo. While the 
current census of the main belt becomes substantially incomplete at a diameter of about 2 km, we 
have detected dozens of new asteroids in the 200-500 meter size range with an 0.9m telescope. 
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